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TECHNICAL MEMORANDUM 1367 


HEAT TRANSFER, DIFFUSION, AND EVAPORATION^ 

By Wilhelm Nusselt 

Although it has long been known that the differential equations of 
the heat-transfer and diffusion processes are identical, application to 
technical problems has only recently been made. In 1916 it was shown 
(ref. l) that the speed of oxidation of the carbon in iron ore depends 
upon the speed with which the oxygen of the combustion air diffuses 
throrigh the core of gas surrounding the carbon surface. The identity 
previously referred to was then used to calculate the amount of oxygen 
diffusing to the carbon surface on the basis of the heat transfer be- 
tween the gas stream and the carbon surface. Then in 1921, H. Thoma 
(ref. 2) reversed that procedure; he used diffusion experiments to de- 

vnr-. ^ Tl_ ^ ^ — 4- 1 T ^ -.V. T ^ 1 

.UJ.J.I.XV., UX Ctixox X J.^LXUO • X\^(.;CL10_L.jr J_A.^XXX XOV^XX CX • J XXCXD 

extended this work by experiment. A technically very Important appli- 
cation of the identity of heat transfer and diffusion is that of the 
cooling tower, since in this case both processes occur simultaneously. 

A relation obtained in the course of such an analysis was given by 
Lewis (ref. 4) and checked by Robinson (ref. 5), Merkel (ref. 6), and 
Wolff (ref. 7). 

In the following it will be shown that more accurate equations 
must be substituted for those used in the previous investigations of 
the relation between the quantity of matter exchanged by diffusion and 
the quantity of heat transferred by conduction. 

A rigid body having a uniform surface temperature T^ is cooled 

by an air stream having a speed wq and a temperature Tq. According 
to Fourier, a quantity of heat 

d^Q » dP dt (1) 

then flows into the air stream in the direction n normal to the sur- 
face element through a surface element dF in the time dt. In this 


^"Warmeubergang, Diffusion und Verdunstung." Z.a.M.M. , Bd. 10, 
Heft 2, Apr. 1930, pp. 105-121. 
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equation, ^ is the instantaneous teii 5 )erature gradient at the surface 

in the direction of the normal, and X is the thermal conductivity of 
air. If equation (l) is applied to an element of space in a gas stream, 
the following differential equation is obtained for the ten 5 )erature 
field in an air stream: 


- -.n2„ 

rcp ^ = XV T 


( 2 ) 


where^ 

Y density of air 

Cp specific heat of unit mass 

X thermal conductivity of air 

Therefore, the Navier-Stokes equations of motion involving u, v, and 
w are to be understood. 

For the diffusion problem, air can be treated approximately as a 
homogeneous body, since the moleculstr weights of its components 
nitrogen and oxygen are only slightly different. Air is now considered 
to be mixed with the diffusing gas, for exeimple, ammonia or water 
vapor. The concentration of the water vapor at an arbitrary point in 
the air stream is designated cj that is, for example, there are c 
kilograms of ammonia in 1 cubic meter of air. In the similarity ex- 
periment of Thoma (ref. 2), the surface of the rigid body was made of 
blotting paper saturated with concentrated phosphoric acid. Amimonia is 
absorbed very actively by such a surface, so that the partial pressure 
and, hence, the concentration of ammonia is very small. In a cooling 
tower the diffusion stream proceeds in an outward direction from the 
water drops. The vapor pressure and, hence, the water -vapor concen- 
tration at the surface of the drop is accordingly dependent upon the 
water temperature. If p^ is the partial pressure of the diffusing 

gas at the surface of the body and pq that in the air stream, then 

the driving force of the diffusion stream is the partial-pressure 
difference Pw ~ Pq* 

^The following nomenclature is used herein in the case of variables 
dependent upon time and the three coordinates x, y, and z; 

® ^ u to ^ , to ^ „ ST 

dt ot ox oy oz 

V% = ^ ^ ^ 

^y2 
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If now, again, n is a running coordinate representing the normal 
to the surface element dF, then the quantity of vapor diffusing through 
\ the surface dF is given by the basic relation of Fick; 

d^Gj^ » -k || dF dt (3) 

in which 

c vapor concentration, kg/cu m 

bc/bn gradient of vapor concentration in direction normal to 
surface 

k diffusion coefficient, sq m/sec 

If an element of space is taken in the air stream, exchange of 
vapor between such an element and its environment occurs partly through 
diffusion and partly through streaming (sensible motion). This fact 
leads to the differential equation: 






^ = kV^c (4) 

dt 

The agreement between equations (l) and (3) as well as between (2) 
and (4) is immediately recognizable. Accompanying a calculation of 
ten?>erature field from diffusion field or vice versa, there must also 
exist, however, an equivalence of boundary conditions. The temperatures 
Tw and Tq of the heat transfer correspond to the gas concentrations 

c-y, and cq of diffusion. But, while during heat transfer the gas 

velocity is zero at the surface of the cooled body, there exists at 
that point, in the case of diffusion, a finite- gas-velocity component 
normal to the body surface. This difference is easy to see if a one- 
dimensional diffusion process is studied in a tube under steady-state 
conditions , Suppose a tube I meters long is filled with air and 
ammonia. Suppose also that by certain experimental means the concen- 
trations at the ends of the tube are held at different, although con- 
stant, values. Then a quantity kilograms per hour of ammonia 

diffuses through the tube in one direction, and a quantity of air G£ 
in the other direction. If x is a running coordinate, Cj^ the 
ammonia, and 03 the air concentration, then 

dCi 

Gi = -ki -jj- 


®2 


dco 

^2-d 


and 


( 5 ) 

(5a) 


4 


NACA TM 1367 


pressures 

of ammonia and of air are Pi and 

Pg, respectively, then 


dPl 

GiRiT = -ki IGT 

(6) 

and 

dpg 



G2R2T = ^2 

(6a) 


Integration of equations (6) and (6a) yields^ 

GiRiT = -kiPi + Cl 


and 


^2^2^ * ^2P2 ^2 


( 7 ) 

(7a) 


hence a linear variation of psurtial pressures. If gravitational in- 
fluence is now ignored, the total pressure pq in the tube is constant. 

Hence, 

Pi + P2 * Po (8) 

It then follows from equations (7) and (7a) that 

ki = k2 = k (9) 

and 

GiRiT = G 2 R 2 T = Vpo (10) 


that is, the same volume, evaluated at the total pressure, diffuses in 
both directions.^ Moreover, with the notation of figure 1, 

V » I (Pi' - Pi") = I (P2" - P2') 


and 



( 12 ) 




0 


3 

These are obviously in error; x lacking. 

^1 and kg are equal, but it does not follow from this argument. 


^9sy 
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as well as 


G2 - I (C2" - C2') (12a) 

Therefore, according to this exan5)le, the partial pressure and the con- 
centration of ammonia vary linearly along the tube axis. Now, the 
experimental situation in the case of the cooling tower and in the 
diffusion research of Thoma (ref. 2) was otherwise. In the case of the 
tube, the boundary conditions are different at one tube end. If that 
end is closed with blotting paper saturated with phosphoric acid, there 
is set up, for the case of linear (axial) diffusion, the equivalent of 
Thoma 's experiment; or, if that end is closed and the bottom covered 
with water, a situation corresponding to that of the cooling tower is 
obtained. Further, Stefan (ref. 8) and Winkelmann (ref. 9) determined 
the diffusion coefficient k for the diffusion of water vapor in air 
with this arrangement. It will now be shown that in this instance 
linear variation of the partial pressure along the tube axis is not 
attained in the case of the stationary diffusion stream. Once again, 
air and ammonia, or water vapor, are counterdif fusing. Since, however, 
one end is impermeable to air, no transport of air along the tube can 
take place. Because of the gradient of partial pressure of air, air 
must diffuse in the axial direction. This molecular air transport must 
work against a convective air transport; that is, a sensible flow of 
gas in the direction of the tube axis must occur. This flow will be 
represented by the symbol u. Then, for the flow of air the equation is 

dco 

Gg *= 0 = ucg + k (13) 


or, after introduction of the partial pressure, 

up2 = - k (l3a) 

For the ammonia stream the corresponding equations are 


(14) 


(14a) 


dcj^ 

= uci - k -33^ 

and 

<iPl 

Gl^lT = upi - k 
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Now again 


Pi + P2 “ PO 

It follows then from equations (l3a) and (l4a) that 


k ^ 

Pq - Pi dx 


( 8 ) 


(15) 


If equation (15) is substituted for u in equation (l4a), the following 
expression is obtained for the partial pressure of ammonia: 


kpo dpi 
GlRlT = - Pq - Pi ^ 


(14b) 


the solution of which, with regard to the boundary condition Pi = Pi' 
at X = 0, is 


Pq “ Pi 
Po - Pi' 


GiRiT 

^Po 


(16) 


Therefore, no linear variation of partial pressure along the tube axis 
occurs in this case, but rather a logarithmic. If pi » pi" at x = 2, 


G 


1 


^^PO , Po - Pi" 

RlTI PO - Pi' 


(17) 


This equation, given earlier by Stefan (ref. 8), shows that the rate of 
gas diffusion through the tube is no longer proportional to the partial- 
pressure gradient. First, if pi is small compared with Pq, it is 
again true that 


c ' - c. 


k Pi' - Pi" n 

i — r 


(18) 


From equations (l5) and (16) it follows that 


u 


RnT 

a,_.vo 


( 19 ) 
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That is, the gas velocity u is identical with the volume of ammonia 
gas streaming through a unit of surface of the tube cross section per 
unit time, as measured at the total pressure Pq. Hydrodynami cally , 

therefore, there exist in the ammonia problems sinks of strength u 
corresponding to evaporation from the source of water. Since, in the 
case of heat transfer, the velocity is at that point zero, an exact 
similarity between heat transfer and diffusion is nonexistent in the 
experiment of Thoma. Only at very small values of hence at low 

partial pressures of ammonia, can u be taken as zero. Only then does 
similarity of the boundary conditions also exist. This will be dis- 
cussed later. Similarity requires, however, identity of the hydro- 
dynamic equations as well. In the experiment of Thoma, the temperature, 
in the case of the diffusion work, is uniform throughout; while, in 
the case of heat transfer, it varies. In this instance, the air 
density varies with temperature. In the case of diffusion, the gas 
densities vary with the concentrations. Similarity obtains, therefore, 
only if 


^w _ T’v 
^0 "'^0 


( 20 ) 


In equation (20) Xv Tq 'ttie densities of the ammonia-air 

mixture at the wall and in the free gas stresim. Further, the variation 
of viscosity with temperature and ammonia concentration must still be 
considered. Since, however, as was shown previously, similarity is 
possible only at small ammonia concentrations , the experiment gives the 
heat -transfer coefficients only at very small temperature differences . 


The similarity conditions are more favorable in the case of the 
cooling tower. In this instance, a hot water surface is cooled by a 
cold dry air stream. Accordingly, heat flows from the water to the 
air, and simultaneously the water that diffuses into the air evaporates. ^ 
In this instance, heat transfer and diffusion are, accordingly, super- 
posed on each other under the same stream conditions. The same hydro- 
dynamic equations and boundary conditions are associated with eqviations 
(2) and (4). 


In the following discussion, the relations that follow from 
similarity considerations in the case of Thoma 's experiment and in the 
case of the cooling tower are derived. 


^Stated precisely as written by Nusselt. 



8 


NACA TM 1367 


I. MDDEL EXPERIMENT OF THOMA 

A body having a surface F and teinperature is cooled by an 

air stream having a speed wq and a temperature Tq* The heat loss of 

the body has the value Q kilocalories per hour. A diffusion experi- 
ment is now carried out by using a body of the same shape. The air 
streeun has the same speed wq and is mixed with ammonia to the concen- 
tration CQ* Through the diffusion process, which occurs at the body 
surface, the ammonia concentration has already attained the value c^. 

The diffusion experiment shows that the surface absorbs G kilograms 
of ammonia per hour. What relation exists between Q and G? This 
question is discussed in the following on the basis of different theories. 


Similarity Theory 

The similarity theory of reference 10 leads to the following ex- 
pression for the rate of heat exchange: 

« > XL{T„ - Hj) f , 

where 

X thermal conductivity of air 
L dimension of body 

f initially unknown function of two dimensionless fractions, 
dependent on shape of body 

X density of air 

q viscosity of air 

Cp specific heat of air 

For diffusion, correspondingly, 

G = kL (cq - c^) f 



( 22 ) 
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From these is obtained the relation sought: 



(23) 


It follows from most of the experimental work that the function f may- 
be represented as a product of two functions, that is, 


f » b 



(24) 


T.TVl£i“V*C 


Ulii. "ta* t.AVAWLXi vr WLA WiiX^ 


11 shows that, for the flow of gases through a cylindrical tube. 


and 


Two\ reasons then occur to set the two exponents equal to each other. 

If it is assumed that the velocity profile over the cross section is 
independent of density, the velocity components u, v, and w are 
proportional, in the differential equation of heat conduction, to the 
stream velocity at the center Wq. In this equation, the fraction 

wycp/x as a factor can be taken out. In the function f it must 
then follow that m = n. Since T^c^g/x varies only slightly among 

the different gases, the influence of the magnitude of n on Q is 
slight. On that account an m « n was chosen, and in that manner a 
very simple kind of equation was obtained. Such a choice was expressly 
limited to gases, because it appeared, since the experimentation of 
Stanton with water had given a greater value for m than Nusselt had 
found for gases, that the power form is valid only in a narrow range 
of values of From the more recent work of Sonnecken and 

Stender (see ref. 12) were obtained, depending on the experimental 
conditions, values of m between 0.72 and 0.91 and values of n be- 
tween 0.35 and 0.50. Merkel (ref. 13) gives for the same research 
m= 0.87 and n = laJz = 0.435. Rice (ref. 14) proposes for flow in 
smooth tubes m « 5/6 a 0.83 and n = l/2 » 0.50. Lately Schiller 
and Burbach (ref. 15) have again grappled with this question. They find 


m a 0 . 786 


n a 0.85 


(25) 
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that the Nusselt relation for gases with m = n is also experimentally 
confirmed for water and, further, support theoretically the equality of 
the exponents. This important question will he discussed later. 

Equations (23) and (24) lead to 



(23a) 


If the heat-transfer coefficient a is calculated with the use of 


then 


Q = (^^(Tw - To) 



(26) 


(27) 


Relation of Thoma-Iohrisch 

In the derivation of his relation whereby the heat-transfer 
coefficient may be calculated from a diffusion experiment, Thoma pro- 
ceeds on the basis that the following condition holds: 


k s 


X 

CpT 


(28) 


In that circumstance the two differential equations (2) and (4) are 
interchangeable . The temperature field is then proportional to the 
concentration field, and it follows that the relation 



^c ^0 ■ *^w 
Sn 


(29) 


holds at any arbitrary point in the fields. From equations (l) and (3) 
the relation sought follows immediately: 


Q ^ X ^w ~ ^0 
G k cq - 


( 30 ) 
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which, using equation (28) , can be also written 



^w - Tq 
•^0 ■ *^w 


(30a) 


If the heat-transfer coefficient a is calculated with equation (26), 
there is obtained 


a = 



(31) 


Equation (28) is not valid in the diffusion experiment with ammonia in 
air, however. 

According to Thoma, it is considered that 


_1_= 1.25 (28a) 

rcpk 

In this case, Thoma uses equation (30a) as well as (3l). As a 
correction factor he multiplies their right sides by the constant of 
equation (28a). He therefore obtains 


Q _ X ~ ^0 

G rCpk ^ Cq - Cw 

from which equation (30) is obtained. Then it follows that 


(32) 


G X 1 
F k cq - c„ 


(32a) 


Equation (27) assumes this form if n =0. 

On the basis of this equation (52a), Thoma derives, from similarity 
and inipulse considerations, the equation for heat transfer: 


Q = gCpTjL(T^ - Tq) (pi 


(33) 
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If this equation is written in the form 


Q = XL(T^ - Tq) ^ gT) 


(33a) 


it is recognized as a special case of equation (21), to the extent that 


f 


fapr ncpg’ \ 

Ven ' \ / 


ncpg 


<p 


X gn 

is valid. From equation (33) it follows that 




which may also he written 


a = 


gn 


X F 6^ 


(34) 


(35) 


(35a) 


Thoma now uses equation (35a) to derive equation (32a) from (31) . The 
author cannot follow Thoma 's reasoning, because if f is replaced in 
equation (23) by Thoma 's relation of equation (34), the following is 
obtained; 


Q = 


G 



(30a) 


and 


a 


G 1 

= F Cq - cw 


(31) 


Thoma 's superaddition leads back, therefore, via equation (34), to 
equation (31) and does not give equation (32a), which was used by Thoma. 

Hitherto it has always been assumed that an alteration of a field 
is effective only in the vicinity of a body. If the whole stream 
cross section is affected, then another definition of a is necessary. 
Set 


2951 


2951 


NACA TM 1367 


13 


dQ = adP(Tw - T) (36) 

wherein T is the temperature in the center of the stream cross section 
that encounters the body sui'face element dF. 

Moreover, the equation 


Q = VoTCp(T2 - Tq) 


(37) 


applies, in which Vq is the air volume swept out by the body in unit 
time, and Tg - Tq the temperature rise of that air caused by the hot 

body. From equations (36) and (37) it follows® that 


Vorc_ T.. - T(^, 
a = — In =— 


F 


^w “ Tg 


(38) 


If 


= € 


then 


^2 - Tq 
Tv - To 

''^orcp 1 

a - p In ^ ^ 


In the diffusion work for the same air volume Vq, there is 
obtained 


(39) 


(38a) 


e = 


' °2 
^0 “ 


(40) 


provided cg is the ammonia concentration behind (downstream of) the 

body. For the first special case, which is presupposed by the validity 
of equation (28), there would be 


C = 6 

and hence the heat-transfer coefficient sought: 


( 41 ) 




6 


It follows only if certain additional assumptions are made. 
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a = !2 I!p in _1_ (42) 

F 1 - e 

In the general case, Thoma and Lohrisch, as in the preceding, follow 
equation (28a) and put 


a = in _i_ (43) 

rcpk r 1-6 

hence, 

a = Z9. ^ In i — (43a) 

F k 1-6 

The correct relation is obtained as in the following . The 
temperature Tg is eliminated between equations (37) and (38), and 


a = !2 I!p in 1 

F Q 

VoTCp(T^ - To) 

Now, according to equation (23a) 

q _ X . 1 

Tw - To " lx J Cq - 


If it is noted that 


G = Vq(cq - cg) 


then 


Q _ X ^0 ' ^2 

Tv - To ^ \ \ J Oq - 0^ 

is obtained, and, therefore, from equations (44) and (40), 


^Q7~^P 


In 


1 



-n 

t 


(44) 


(23a) 


(45) 


(46) 


a 


F 


1 


(47) 
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This equation should therefore replace equation (43a) of Thoma- 
Lohrisch. According to Thoma's research, e is in the vicinity of 0.25. 
If the numerical value, of equation (28a) is inserted in equation (47) 
and, further, the values n = m= 0.6, equation (47) yields a value of 
a some 11 percent smaller than Thoma's equation (43a). 

An additional inaccuracy arises also in the calculation of e 
according to equation (40). Thoma, in that equation, put Cy = 0, 
since he assumed that in consequence of the strong absorption of am- 
monia by phosphoric acid the partial pressure of ammonia in that rer 
gion is zero. Since, however, the ammonia must diffuse through the 
boundary layer, a finite vapor pressure of ammonia must exist at the 
surface. It can naturally be quite small, but it is necessary, first, 
to measure it once. Lohrisch also employed Thoma's experimental tech- 
nique in the case of water vapor, to the extent that he saturated with 
water the blotting paper comprising the body surface. In the calcu- 
lation of 6, he assumed that the vapor pressure at the body surface 
corresponds to the water temperature. In section II, concerning 
evaporation in a cooling tower, it is shown that the vapor pressure 
at a body surface is smaller than the saturation pressure. 


Boundary-Layer Theory 

In its most primitive and simplest form, this theory supposes that 
an air stream flowing past a body may be considered to consist of 
two contiguous but sharply demarcated portions, namely, the boundary 
layer adjacent to the body surface, and the balance of the air stream . 
The one is associated with a laminar flow, the other with a turbulent, 
which in the first approximation is treated as a potential flow. Then 
it is assumed that in the latter a complete equalization of temperature 
or of ammonia concentration occurs. Figure 2 exhibits this distri- 
bution of temperature T and concentration c in the case of heat 
transfer or of a similarity experiment. The width of the boundary 
layer is indicated by y. For heat transfer, then. 


>^(Tw - To) 

y 


(48) 


and for diffusion 

kF(cQ - c.„) 


y 


( 49 ) 
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Accordingly, there is obtained on the basis of elementary boundary- 
layer theory 

i i ~ *^0 (50) 

G " k CQ - 

and also 

a = £ ^ i (51) 

F k cq - c^ 

which is identical with the equation (32b) of Thoma and with equation 
(27) when n = 0. 


Impulse Theory 


Impulse theory is a carry-over of gas -kinetic considerations into the 
domain of turbulent motions of a fluid. In its simplest form it is as- 
sumed, as an explanation of heat transfer, that a volume of gas V having 
the temperature Tq moves from the turbulent fluid stream to the wall, 
where it is heated to a temperature and then brought back to the core 

of the fluid stream. It therefore takes from the wall the heat 

Q = VrCp(T^ - Tq) (52) 

which is given up to the fluid. It nrust at the same time he true that 


VTCp 

F 


(53) 


The volume V contains in the case of diffusion Cq kilograms of 

ammonia. On impact against the body surface, ammonia is absorbed until 
the concentration is c^. Hence, 

G = V(CQ - c^) (54) 

From equations (52) and (54) it follows that 



rc 


'^w - ^0 

P =0 - S/ 


(55) 
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and, 


accordingly 


7 




a 




which follows also from equation (27) when n 


1 . 


(56) 


Comparison of Boundary-Layer and Impulse' Theories 

It will now he assumed that in the turbulent core neither perfect 
equalization of speed nor, therefore, of temperature or concentration 
is attained. On the core side of the boundary layer, the temperature 
Tg and concentration Cg are then, respectively, different from 

Tq and Cq. Within the boundary layer, exchange occurs in accordance 

with boundary -layer theory. In the free gas stream, impulse theory 
applies. 

Then, in the boundary layer. 


Q = I F(Tw - Te) 

and 

G = I F(ce - c^) 

Hence, 

a 

G k \cg - 

In the turbulent stream, on the other hand. 


Q = VrCp(Te - To) 


and 


G = V(cq - C 0 ) 

therefore. 


Q *^e ~ '^0 

G “ CO - Cg 

"^Equation (56) is identical with (3l). 


(57) 

(58) 

(59) 

( 60 ) 
• ( 61 ) 


( 62 ) 
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It can now be assumed that, approximately, 


'^w “ ^e 
Tw - To = 


that is. 


is equal to a constant. 


8 


Also, 


b 

” ^w 

From equations (63) and ^64) it follows that 


T, 


T, 


0 


Tw - To 


= 1 


- a 


and 



CO - <Hr 


In that case, equations (59) and (62) become 


Q _ X a Ty ~ Tq 
G ~ k b cq - 


and 


Q 1 - a Ty ~ Tq 

I = 


D 

Recently, Prandtl (ref. 16) has given for this 


gCpil ji 

en 

X " \ 

wo dr 

1 . f - 1 

)e 

V X 

/ yj WQdr 


in which the value of the parameter e is uncertain 
1.1 and 1.75. 


(63) 

(64) 
(63a) 
(64a) 

(59a) 

(62a) 

quantity the value 
It lies between 
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From the last equations, 






CO 

CO 

oi 


a 


a 


(l-a) 


rcpk 

X 


There is obtained, therefore, the desired relation; 


(65) 


G " k 


a + 


(l-a) 


X ^ 


( 66 ) 


If this is compared with equation (23a), it is seen that the two 
functions can be distinguished only with respect to the dependency on 
the fraction ycpk/x. Since n > 0 and 0 < a< 1, both relations lead 
to increasing values [of Q/G] with increasing values of ycpk/x. In 

the case of gases, only a small range of variation of this fraction, in 
the vicinity of unity, occurs and is, therefore, of significance: If 

agreement of the two relations is demanded, there is obtained as the 
connection between two constants the expression^ 


n = 1 - a 


(67) 


Then, if n = 0.4, a = 0.6. 


II. EVAPORATION OF WATER 

Stefans should be credited with having first recognized that the 
evaporation of water is a problem of diffusion. At the same time, he 
developed the theory of diffusion. It is necessary to distinguish 
among several different cases in connection with diffusion. 

Consider first a quiet surface of water at the same temperature as 
the overlying air. If the relative humidity of the air is less than 
100 percent, water evaporates; that is, superheated water vapor diffuses 
into the air from the water surface. Since, under the same conditions, 
water is lighter than air, an air-streaming occurs. Above the water a 
rising current of air develops that sucks dry air over the water sur- 
face. Through that mechanism the evaporation is increased. If a wind 
blows over the water surface, a further increase in evaporation occ\rrs 
as a consequence of turbulence. If the water temperature differs from 


\his, from a + (l-a) (l+A) = (l+A)*^, 


or a + l + A- a - aA = 1+ nA+... etc . 
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the air temperature, an intrinsic influencing of the evaporation occurs, 
while because of the resulting heat exchange the convection is 
influenced. 


Evaporation in Still Air and in Uniform Temperature Field 

Above a surface of water of area F having a representative 
dimension L, a layer of air exists, the density of which, at some 
distance away, is xo 'tbe specific humidity of which is Cq. At 

the water surface the humidity content of the air is c^. At any point 

above the [water] plane the vapor content of the air increases as a 
consequence of evaporation; it becomes lighter and suffers a lift Z 
in the amount 


Z = 



( 68 ) 


where |i is the apparent molecular weight of dry air and ^bat of 

damp air. For the diffusion field, diffusion equation (4) then applies. 
The Navier-Stokes equations of motion are required for the determina- 
tion of the velocity components appearing therein (eq. 4); in these 
equations, the lift Z appears as an external force acting in the 
direction of negative gravity; and, therefore. 


lo ^ - 7 

g dt “ ^ 




(69) 


In this expression, the air density Xq can be assumed constant. 

A similarity consideration leads to the expression for the concentration 
gradient at the water surface: 


12 . 

dz 




$ 


3/iL 


- l) (Cw 


Co) 




kr 


(70) 


The vapor mass evaporating per unit time from the water surface 
G is calculated according to equation (14). If the value of u accord- 
ing to equation (l9) and the following value of c^^ 

S? - RT 


(71) 
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are substituted in equation (14), then 

_ -IF dc 

Pw dz (^2) 

1 - — 


where p is the total pressure. 


If the concentration gradient is now replaced by its equivalent 
according to equation (70) , the rate of vapor evaporation becomes 



3 Ji_ 

L rin - 1) (cy - cq) 


M 

kr 


(73) 


Evaporation in Wind and in Uniform Temperature Field 

Over a water surface a wind passes whose speed at some distance 
from water has the uniform value Wq. The air and water tei^peratures 
are identical. Gravitational influences can accordingly be ignored if 
the airspeed exceeds several meters per second. Therefore, the ob- 
servations and formulas of section I apply (ref. 17). It follows that, 
by equation (22), the rate of water evaporation is 

G= kL(<=^- Cj) f , m) (32) 

If the assumption is now made that the velocity u arising from 
the evaporation normal to the water surface can be ignored, the function 
f of equation (22) can be taken over from the corresponding heat- 
transfer problem in accordance with equation (21) . The results of 
reference 17 should be considered here. A copper plate heated elec- 
trically to 50° C, and having a dimension on a side of 0.5 meter is 
cooled by an air stream having a temperature of 20° C. For w > 5 
meters per second, the following value was obtained: 

O 7A 

a= 6.14 wq ’ kcal/ (m^ ) (hr ) (°C ) (74) 

It follows that in equation (24), if m * n, 

. .0.78 
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In equation (22) there is obtained, accordingly, 

0.78 

/wqLN 

f = 0.065 J 


(75a) 


Therefore, the expression for the amount of water evaporating from 
a water area of F square meters in an hour becomes 

. . 0.22 

G = 39F (cw - Cq) kg/hr (76) 

where k is the diffusion coefficient according to the research of 
Mache (ref. 18): 


k = 


0.078 


/ T / 

j '1 ^1“ 


in which the total pressure p is to be used. 


(77) 


Evaporation with Heat Exchange in Still Air 

It will now be assumed that the temperature of the water sur- 

face is different from the temperature Tq of the air. Hence, heat 

transfer occurs in addition to diffusion. In this instance, both 
processes are coupled through the resulting air stream to the extent 
that air-lifting is caused by both the lesser specific gravity of the 
water vapor and the heating of the air. Instead of equation (68) of 
section II, the lift Z per cubic meter has the value 

Z = ror(T - Tq) + (c - (68a) 

Here T and c are the temperature and specific hiomidity at any place 
in the field, y"q is the specific gravity of dry air at a great dis- 
tance from the water surface at which the conditions pq and Tq pre- 
vail, and r is the coefficient of expansion of air: 

r = i/t (78) 

In the case of excesses of temperature that are not too great. 


r = 1 /Tq 


(78a) 
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approximately. Then 


Z = ^(T - To) + (c - cq)(^ - l) (68h) 

becomes the value of the body force in the equation of motion (69). 
Moreover, equation (2) of heat conduction and equation (4) of diffusion 
must be used. With the abbreviations 

, - To) 

2 

en^To 

c = 

gCpT) 

„ ik ■ 0 <=» - s 

^ 2 
rrrt ^ 




\ (79) 


and 





similarity considerations lead to the relations 


G = kL(c^ - Cq) <^1 (B,E,C,D) 


(80) 


and 


Q = XL(Tw - Tq) ^2 (B,E,C,D,) (80a) 

where $ and #2 initially unknown functions of the variables 

B, E, C, and D. 


The formulation is significantly simpler and clearer if it is 
assumed that equation (28) applies, that is, 



for then 


C = D 


(28) 

(28a) 
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In this C 6 LS 0 T &nd. c 8 . 1*6 pi*oportion&l to 68-ch othor* Ecj^n&tion (29) 
applies and, then. 


Sl 

G 


Tv -^0 


X(T^ - Tq) 

" X(Cv - - 


(81) 


The exchanged heat and evaporated moisture stands therefore, in a very 
simple relation to each other. If the values of G and Q according 
to equations (80) and (80a) are substituted in equation (81), 


$ 


1 



(81a) 


which is valid only if equation (28) obtains. For the determination of 
the function # there is introduced, as a matter of expediency, a new 
dependent variable; 

a = T - To + ^ - 1 ) (c - Co) (83) 

which then leads to the differential equations 


ro dw ro 2 

g dt - To ^ 


and 


rocp § = xv2s 


(83) 


from which, through consideration of similarity involving the gradient 
at the water surface, it follows that 


• 5E = 


J-w 




ik ■ 0 





(84) 


This gradient can also be calculated by the use of equation (82), from 
which it follows that 


as _ St , Tq /p _ \ Sc 

Sn Sn Tq \^i ) ^ 


( 85 ) 
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If the divergence ("speed of expansion") u at the water surface is 
considered negligible, equations (l) and (3) apply: 

and 

a = -kF If (3) 


From equations (l), (3), (81), (84), and (85): 


and 


G «= f ^b+E), dJ 

Q = y 1 Q^ f j(B+E), cj 


( 86 ) 


The second of the equations (86) must also apply for the case 
•^w “ *^0 " therefore when E = 0. The equation then goes over 

into the usual form for heat transfer. If the latter is known, the 
rate of evaporation G can thus be calculated. 


On the basis of heat-transfer research (ref. 19), for large values 
of B/C, the following can be written: 


Q = 



- Tq) - CiX 




(T^ - Tq)^"^^/^ 


(87) 


in which the coefficient C-j^ is dependent upon the form and orientation 
with respect to gravity. 

With equation (87), one now obtains from equation (86) 



( 88 ) 
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It should be observed, above all, that this equation is valid only 
if equation (28) is valid, and, hence, C = D. The values k, X, y, qj 
and Cp pertain to the vapor-air mixture and, indeed, are to be taken 
as mean values over the whole field. The term k is given by equation 
(77). The thermal conductivity of air \ is in the following form in 
reference 10: 


X = 


0.00167(1 + 0.000194T) Vt* 


1 + 


117 


kcal/ (m) (hr) (°C) 


(89) 


Since the thermal conductivity of water vapor in the germane tem- 
perature range is only slightly less than that of air, the thermal 
conductivity of the vapor-air mixture can be assumed equal to that of 
air. The terms y ^•iid Cp are to be calculated according to the 

relation appropriate for a mixture of gases. Thus, there are obtained 


at t = 0 


X 

krcp 


0.87 



(90) 


This fraction is different from unity for this vapor-air mixture. 
Accordingly, equations (80) and (80a) apply. Since equations (80) and 
(80a) become equations (87) and (88) in the limiting case C = D, in 
general, in the first approximation, the following can be written: 


and 


4 



(91) 


(92) 


Then the following is obtained: 


G _ k 
Q " X 



(93) 


G _ k a| ^ - Co 

Q X y kCplo - To 


( 94 ) 


or 
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CO 

CO 

cn 


If the temperature within the domain falls below the dewpoint, fog 
formation is initiated and the equations become invalid. 

Equation (87), along with equations (88), (91), and (93), is 
applicable only at large differences of temperature and water-vapor 
concentration. For small differences, the function f(B/c) can be 
expressed only by means of a graphical representation. By such a 
representation, it can be shown that for B + E = 0, f approaches a 
constant value. 


If the air temperature Tq is greater than the temperature of the 
water surface it is possible, for finite values of Tq - and 

Cw “ cq, that B + E = 0. In this case, no convective streaming occurs, 
but rather only a molecular transport of heat and vapor. Hence, for 


B = -E 

or 

•^0 - ^w ^ ^ ^ 

°w - <=0 ^0 V^^l ' 

there is obtained 

G = C2kL(cw - Cq) 

and 

Q = C2\L(Tq - T^) 

Therefore , 


G _ k Cw - Cq 
Q X Tq - T^ 


(95) 
(95a) 

(96) 

(97) 

(98) 


The coefficient C2 is dependent upon the shape of the water surface. 


Evaporation with Heat Exchange and Air Flow 

As in the case of Evaporation in Wind and in Uniform Temperature 
Field, a wind having a speed Wq flows by a water surface. However, 

the water temperature and the air temperature Tq are now dif- 

ferent. With a small partial pressure of water vapor assumed, equations 
(21), (22), and (23) of section I apply here as well. Since, however. 
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at 50° C the partial pressure of water vapor has already attained a 
value of 0.125 atmosphere, they must he modified. For the Reynolds 
number the following will, for brevity, be used; 


Re 


Lworo 

ST 


(79a) 


On the basis of similarity theory one has again, first of all, the 
relations 


and 


dT 

■ SH = 


Tw - To 
L 


f (Re,C) 


dc 

■ = 



f (Re,D) 


(99) 


(99a) 


In these, in the case of a smooth water surface such as a tank,^® L is 
the principal dimension; and for a water drop, it is the diameter. At 
the water surface there exists, further, between the rate of evaporation 
and the concentration gradient the following relation; 


G 

F 


V 

ucv - k 3^ 


(14^ 


where, according to equation (19), 


u 


GRiT^ 

Fp 


(19a) 


and Rj_, in this equation, is the gas constant of water vapor, and p is 
the total pressure. There is then obtained 


uc^ 


^Pw 

fT 


( 100 ) 


in which p^ is the partial pressure of water vapor at the water surface. 
If the value of equation (lOO) is used in equation (l4a), the concen- 
tration gradient at the water surface becomes 


do _ g (l 
3n kF 


( 101 ) 


10 


Or pond or pool. 
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By the same reasoning, the following applies for heat transport: 


f - '^TCp (Ty - Tq) - ^ ^ (I4t) 

and, with equation (19a), 

• i [« ■ ° S;; =p<’'w - ’’o)] (“2) 

From equations (99), (99a), (lOl), and (102), the relation sought 
between Q and G is 


Q / T ^^w \ “ To)f(Re,C) a , . , . 

G " \ ^ ■ P / h(c^ - CQ)f(Re,D) - Tq) (103) 

If the value of f given by equation (24) is substituted in 
equation (103), 

8- (~.i 

From equations (lOl) and (99a) there is obtained, with equation 


(24), 


G = 


bkF(c^ - cq) pg 


m 


- t) 

and, at the same time, from equation (I03a), 


ID^ 


(104) 


Q = 


bF(T^ - Tq) 






“1 f-F)d 


Re 


m 


( 105 ) 


If Cy = Cq, no evaporation occurs and equation (105) becomes 
equation (24). 
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Vapor Pressure at Water Surface 

It is natural to assume that, at the evaporating water surface, the 
vapor pressure is equal to the pressure at the saturation value corre- 
sponding to the water temperature and that, therefore, P-y^ = Ps* 

At the same time, the vapor concentration at the water surface then 
becomes 


Cw ^ (106) 

that is, it is equal to the saturation density. That this, however, is 
not the case has already been conjectured by Winkelmann (ref. 9) and 
then demonstrated by Mache (ref. 18), who found, on the basis of a 
thorough research on evaporation in a cylindrical tube, that the follow- 
ing relation exists between the rate of evaporation and the pressure in 
question: 


|=Kb(Ps-Pw) (107) 

that is, the vapor pressure over the water surface, during evaporation, 
is always smaller than the saturation pressure corresponding to the 
temperature of the water surface. The coefficient Kq is a temperature 

function that unfortunately has not yet been precisely determined. If 
the density of vapor instead of the partial pressure is introduced in 
equation (l07). 


|=Pl(c"-Cw) ( 108 ) 

where c" is the saturation concentration of water vapor at the water 
temperature Tw> and is a constant dependent upon temperature, 

which, according to the researches of Mache, assumes values dependent 
upon the water temperature as indicated in the following table of values: 


tw, °c 

Kq, 1/hr 

3l, m/hr 

92.4 

0.0086 

148 

87.8 

.0080 

133 

82.1 

.0084 

108 

27.5 

.084 

925 


Unfortunately, it is precisely in the technically important temperature 
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range between 0° and 50° C that only a single experimental value is 
available. 


Application to Psychrometer of August 

For many technical applications of the diffusion relations pre- 
viously developed, it is appropriate to calculate the evaporation 
coefficient on the basis of the heat-transfer analog, as has already 
been done in the treatment of burning and vaporization of the carbon in 
iron ore (ref. 1). Set 


G = 3F(c" - cq) kg/J^ (109) 

from which the dimensions of the evaporation coefficient are 


3 = m/hr (110) 

In the case of August's psychrometer, evaporation takes place from 
a moist thernometer in still air. In such an instance, equation (91) 
applies in the calculation of the mass of water evaporating per unit 
time. If, further, an evaporation coefficient Pg is inserted, where 

( 111 ) 

then, equation (91) becomes 

G = P2F(<V - (91®^) 

Further, equation (108) still applies. If the unknown concentration 
c^ is eliminated between equations (108) and (91a) , the evaporation 
coefficient p in equation (109) becomes 




P h h 


( 112 ) 


When the heat transferred along the stem of the thermometer is ignored, 
the heat balance of the wetted thermometer may be expressed as 


Q * G 



(113) 


In this expression, Tw is the temperature of the wetted thermometer, 
Tg the saturation temperature corresponding to the partial pressure 


32 


NACA TM 1367 


of the water vapor at the surface of the wetted thermometer, r the 
heat of vaporization at the pressure p^, and Cp the specific heat of 
the limiting curve [?] at the same pressure. 

If Tq is the environment temperature, as measured with a dry 
thermometer, the heat extracted from the surroundings is 


Q = oF(To - Tt^) 


(114) 


The heat-transfer coefficient a consists of two parts, several terms 
of a sum covering the heat conveyed by thermal conduction to the 

thermometer, and a portion as that gives the magnitude of heat radia- 
tion. The latter is 


as = 



where Cg 


is the radiation coefficient of water. 


that is, 


Cg * 3.35 kcal/(m2)(°C^) (116) 

and Tq' is the mean temperature of the fixed body surrounding the 

wetted thermometer with which radiation is exchanged. It is certainly 
approximately equal to the ambient -air teinperature Tq, yet surely not 

quite precisely equal. Herein, under certain conditions, exists a not 
unimportant source of error in psychrometry. This source of error can 
successfully be eliminated (as was communicated to the author by 
Dipl.-Ing. Kaissling) by surrounding the wet thermometer by a radiation 
shield, which consists, as does the wet thermometer itself, of a 
wetted surface. 

If Tq' = Tq and the attainment of room temperature is assumed, 
equation (115) becomes, approximately,-^^ 

ttg = Cg (115a) 


^^Dimension time"^ apparently missing in equation (116). 

12 

This does not seem to be correct* 
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The heat -transfer coefficient is dependent upon the flow conditions 

in the vicinity of the thermometer. If the psychrometer is hanging in 
a region in which the air is quiet, then, according to equation (92), 


Ob 



B + E 
G 


(117) 


in which the constant is dependent upon the configuration of the 

thermometer well. For a cylindrical well of height H, 


Cl = 0.83 


(118) 


In equation (116 [117]), L is replaced by H. 


If the wetted thermometer is placed in a current of air, there is 
obtained, for example, for a plate-shaped thermometer, the following 
relation (ref. 20): 


X 

= 0.069 L Re + 0.83 




B + E 


-0. 00002 8Re 


(119) 


wherein it is supposed that the wind flows along the thermometer well 
in a horizontal direction. 


If, in the energy equation (113), the value of G from equation 
(109) is inserted, and that of Q from (114), the following psy- 
chrometer formula is obtained: 

j^r + Cp(T^ - Ts)] (c" - Cq) = (Ob + Cg) (Tq - T^) (120) 


For the diffusion constant of a plate-shaped, wetted thermometer, 
following is obtained with equations (119), (lll)> and (117): 


8g = 0.069 


lit) 


0.78 


+ 0.85 j 




-0.000028 


B + E 


it) 


0.78 


the 


( 121 ) 


Equation (120) gives a decrease of the psychrometer constant with 
increase of airspeed, which has been well confirmed by the investi- 
gations of Edelmann, Sworykin, and Recknagel. 
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Application to Theory of Cooling Tower 

In this [ apparatus ] , finely divided warm water trickles downward 
and is cooled by a rising current of cold air. If, at some point, wq 

is the relative speed of the water and air with respect to each other, 
equations (103), (104), and (105) apply for the heat transfer and 
evaporation. If the heat-transfer coefficient a and the evaporation 
coefficient Pg now evaluated according to equations (26) and 

(92a), 

and, with equation (112), 



^fe,terial on the technical applications of the formulas here pre- 
sented will soon be published elsewhere. 


Translated by H. H. Lowell 
National Advisory Committee 
for Aeronautics 


1 -Z 

Since the transmission of the original manuscript to the editor' 
office on May 16, 1929, the following papers have appeared: E. Schmidt 

Verdunstung und Whrmeubergang, Gesundheitsing. , 1929, p. 525.} R. 
Mollier, Das ix-Diagramm fur Dampfluftgemische, Stodolafestschr. , 
Zurich, 1929, p. 438} H. Thiesenhusen, Untersuchungen uber die 

Wasserverdunstungsgeschwindigkeit in Abh'dngigkeit von der Temperatur 
des Wassers, der Luftfeuchtigkeit und Windgeschwindigkeit, 
Gesundheitsing., 1950, p. 113. 
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